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Summary.  The  Mark  III  interferometer  is  an  operational  long 
baseline  stellar  interferometer  on  Mt.  Wilson  with  four  possible 
baseline  configurations  from  9  m  NE-SW  to  20  m  N-S.  The 
interferometer  was  designed  to  be  a  highly  automated  astronom¬ 
ical  instrument  to  measure  stellar  positions  and  diameters  to  a 
magnitude  limit  of  seven.  Initial  fringe  observations  were  made 
in  September  1986  with  a  12-m  N-S  baseline.  In  the  following 
months,  semi-automated  astrometric  and  stellar  diameter  mea¬ 
surements  were  also  made.  This  paper  describes  the  hardware  and 
software  components  of  the  instrument  and  its  operational 
characteristics. 

The  interferometer  has  several  novel  features.  One  is  the  use  of 
optimal  estimation  and  control  algorithms  (e.g.  Kalman  filters)  in 
the  control  loops.  Another  is  the  ability  to  operate  both  as  a 
closed-loop  phased  interferometer  and  eventually  as  an  open-loop 
or  absolute  coherent  interferometer.  High  thermal  stability  and 
mechanical  accuracy  should  permit  the  instrument  to  point  blind 
at  an  astronomical  object  and  maintain  optical  path  equality  to 
within  the  limits  set  by  the  atmosphere.  In  this  absolute  inter¬ 
ferometric  mode  of  operation,  it  should  be  possible  to  observe 
faint  astronomical  objects  that  are  too  dim  for  phase  tracking.  In 
theory,  measurements  of  amplitude,  group  delay,  and  closure 
phase  will  be  possible  to  14  mag. 
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I.  Introduction 

The  Mark  III  interferometer  is  the  third  of  a  series  of  stellar 
interferometers  at  Mt.  Wilson  Observatory.  The  first  instrument 
demonstrated  white-light  fringe  tracking  in  1979  (Shao  and 
Staelin,  1980).  The  second  instrument  observed  fringes  in  1982 
and  was  used  as  a  technology  test  bed  for  astrometric  measure¬ 
ments  until  1984  (Shao  et  al.,  1987).  Several  devices  critical  to  the 
operation  of  a  stellar  interferometer  were  developed  for  the  Mark 
II  interferometer.  One  such  device  was  a  high  speed,  ultra  high 
accuracy,  laser  controlled  optical  delay  line.  An  operational 
procedure  for  rapid  switching  between  stars  was  also  developed. 
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Traditionally,  the  most  difficult  task  of  a  stellar  interferometer  has 
been  finding  and  tracking  the  central  fringe.  In  rapid  switching, 
this  procedure  is  totally  automated.  Thus,  operation  of  the  Mark 
II  helped  identify  the  technical  problems  of  making  astrometric 
measurements  with  a  long  baseline  interferometer. 

Like  the  previous  instrument,  the  Mark  III  interferometer  was 
also  designed  to  be  a  test  bed  for  new  techniques  and  technologies 
relevant  to  optical  interferometry  and  synthetic  optical  aperture 
synthesis.  However,  the  principal  purpose  of  the  Mark  III 
interferometer  was  to  initiate  a  test  program  of  fundamental 
astrometric  measurements  to  demonstrate  a  significant  improve¬ 
ment  in  accuracy.  Another  goal  was  to  initiate  a  program  of 
accurate  stellar  diameter  measurements.  Most  important  of  all, 
the  instrument  was  designed  to  be  astronomically  productive. 

Long  baseline  stellar  interferometers  have  traditionally  been 
considered  one  of  the  most  technologically  challenging  astronom¬ 
ical  instruments.  In  the  Mark  III  interferometer,  the  goal  has  been 
to  build  an  instrument  that  is  reliable,  easy  to  operate,  and  capable 
of  extremely  accurate  astronomical  measurements.  In  order  to 
achieve  these  goals,  a  number  of  active  subsystems  were  in¬ 
corporated  into  the  design  of  the  instrument.  In  implementing 
these  subsystems,  the  design  philosophy  involved  keeping  the 
mechanical  and  optical  systems  as  simple  as  possible  and  shifting 
the  bulk  of  the  complexity  into  software. 

Although  the  Mark  III  was  originally  designed  as  a  closed  loop 
phased  interferometer,  operational  experience  has  shown  that  it 
has  sufficient  accuracy  and  stability  to  operate  as  an  absolute 
interferometer.  An  absolute  interferometer  is  one  whose  mechan¬ 
ical  accuracy  and  stability  in  open  loop  operation,  without  fringe 
tracking,  is  sufficiently  high  that  optical  path  errors  are  dominated 
by  atmospheric  turbulence.  This  level  of  accuracy  and  stability, 
usually  associated  with  the  primary  mirror  of  a  large  astronomical 
telescope,  is  the  level  of  accuracy  exhibited  by  the  Mark  III 
interferometer  with  the  12-m  baseline.  Operation  with  the  20-m 
baseline  is  expected  to  exhibit  the  same  level  of  accuracy  and 
stability. 

The  ability  to  operate  as  an  absolute  interferometer  makes 
possible  three  new  modes  of  interferometric  observations  for  long 
baseline  interferometers.  They  are  photon  starved  amplitude 
measurements,  photon  starved  group  delay  astrometry,  and 
photon  starved  closure  phase  measurements.  Photon  starved 
amplitude  measurements  have  long  been  a  part  of  speckle 
interferometry  with  single  apertures.  The  importance  of  photon 
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Fig.  1.  The  Mark  III  interferometer  on  Mt.  Wilson 


starved  observations  is  that  they  make  possible  astrometry  and 
aperture  synthesis  of  faint  extragalactic  objects.  These  three  types 
of  measurements  are  discussed  after  the  description  of  the 
instrument. 


2.  The  instrument 

The  interferometer  is  located  in  California  on  Mt.  Wilson  ap¬ 
proximately  80  m  east  of  the  100-inch  telescope,  near  the  site  of  the 
Mark  II  interferometer.  Figure  1  shows  a  view  of  the  site.  Six  and 
ten  meters  north,  south,  and  east  of  the  central  building  are  the 
siderostat  piers,  which  are  each  constructed  of  steel-reinforced 
concrete  with  a  4  x  4  foot  base  and  a  height  of  8  feet,  4  feet  of  which 
is  above  ground.  Siderostats  are  currently  mounted  on  the 
innermost  piers.  Huts  provide  weather  protection  for  the  sidero¬ 
stats.  These  have  hinged  roofs  which  swing  back  exposing  the 
siderostats  to  the  sky.  Currently,  there  are  vacuum  pipes  between 
the  huts  and  the  central  building  to  eliminate  local  seeing  effects 
caused  by  heat  in  the  main  building.  Initial  operation  involved 
only  the  12-m  baseline  connecting  the  north  and  south  siderostats; 
however,  the  east  siderostat  can  be  used  in  place  of  the  north 
siderostat  for  a  9-m  NE-SW  baseline.  The  fourth  window  on  the 
central  building  is  for  future  expansion.  In  1988,  we  plan  to 
operate  on  the  outer  set  of  piers  shown  in  Fig.  1,  with  20-m  N-S 
and  16-m  NE-SW  baselines. 

Since  the  instrument’s  principal  purpose  is  large  angle  astro¬ 
metry,  attention  was  paid  to  construction  practices  that  affect 
astrometric  accuracy.  The  central  building  is  air  conditioned 
during  the  day  so  that  at  night  there  will  be  a  minimum 
temperature  difference  between  inside  and  outside.  The  piers  in 
the  huts  that  hold  the  siderostats  go  down  to  bedrock  and  are  not 
directly  connected  to  the  concrete  pads  on  which  the  huts  are  built. 
The  piers  are  also  insulated  by  two  inches  of  foam.  The  optics 
inside  the  main  building  that  affect  astrometric  measurements  are 
all  located  on  a  super-invar  breadboard,  except  the  mirrors  in  the 
optical  delay  line.  However,  the  position  of  the  optical  delay  line  is 
monitored  by  a  laser  interferometer  referenced  to  the  same  super¬ 
invar  breadboard.  The  laser  is  frequency  stabilized  with  better 
than  10  “  9  long-term  stability. 

The  instrument  can  be  divided  into  five  subsystems:  1)  the  star 
tracker  or  autoguider,  2)  the  optical  delay  line,  3)  the  stellar  fringe 


acquisition  and  tracking  system,  4)  the  laser  metrology  system, 
and  5)  the  siderostat  pointing  and  observation  sequencing  system. 
A  description  of  the  optical  layout  of  the  instrument  is  followed  by 
a  more  detailed  explanation  of  the  five  subsystems. 

2.1.  Optical  layout 

Figure  2  is  an  optical  schematic  of  the  interferometer.  Light  from 
the  two  siderostats  is  directed  toward  the  central  building  by  fixed 
mirrors.  These  fixed  mirrors  are  necessary  to  keep  the  polarization 
vectors  from  the  two  interfering  beams  parallel  for  observation  of 
stars  located  over  all  parts  of  the  sky.  After  the  light  enters  the 
main  building,  the  beams  are  directed  by  piezoelectrically  con¬ 
trolled  mirrors  towards  the  vacuum  delay  lines.  After  the  light 
reflects  from  the  delay  line  retroreflectors,  the  beams  are  com¬ 
bined  at  the  beamsplitter  B.  After  the  beamsplitter,  part  of  the 
light  is  directed  into  optical  fibers  which  feed  GaAs  phototubes 
located  some  distance  from  the  central  breadboard,  and  part  of 
the  light  is  directed  toward  the  star  tracker.  Other  optical 
components  include  the  laser  interferometers,  and  alignment 
lasers  and  telescopes. 

2.2.  Star  tracker 

The  star  tracker  has  been  described  in  detail  previously  (Clark  et 
al.,  1986;  Clark,  1986)  and  we  will  provide  only  a  brief  summary 
here.  The  star  tracker  system  is  designed  to  keep  the  two 
interfering  wavefronts  parallel  to  within  a  small  fraction  of  an 
arcsecond.  As  with  all  servo  systems,  the  star  tracker  is  composed 
of  an  error  sensor,  a  controller,  and  an  actuator.  For  the  Mark  III 
interferometer,  two  different  error  sensors  were  built  because  a 
key  component  (a  photon  counting  camera)  for  the  original  design 
was  not  delivered  before  initial  observations  were  made  in 
September  1986.  Only  the  star  tracker  design  that  will  be 
permanently  attached  to  the  interferometer  is  described. 

An  optical  schematic  of  the  star  tracker  is  shown  in  Fig.  3. 
Light  from  the  star  in  a  3-inch  beam,  intercepted  by  the  two 
siderostats,  is  combined  at  a  beamsplitter.  Light  from  the  south 
siderostat  passes  through  an  annular  wedge  that  deflects  the  outer 
portion  of  the  beam  by  45";  the  inner  2-inch  part  of  the  beam  is 
undeflected.  After  the  light  is  combined  at  the  beamsplitter,  an 
annular  mirror  at  one  output  of  the  beamsplitter  directs  the  light 
in  the  annulus  towards  a  telescope  and  imaging  detector.  The  light 
in  the  central  2-inch  core  from  both  sides  of  the  beamsplitter  is  fed 
to  photomultipliers  that  detect  the  fringes. 

When  the  two  components  of  the  light  beam  that  enter  the 
fringe  detection  phototubes  are  parallel  (as  they  must  be  for 
proper  operation  of  the  interferometer),  the  two  components  of 
the  annular  beam  entering  the  telescope  will  be  tilted  by  45" 
relative  to  each  other  because  of  the  45"  wedge  in  the  south  arm  of 
the  interferometer.  Thus  two  star  images  will  appear  on  the  face  of 
the  photon  counting  imaging  detector.  The  optical  arrangement  of 
this  star  tracker  is  slightly  different  from  that  used  for  some  other 
interferometers.  The  present  design  using  one  telescope  for  both 
star  trackers  eliminates  the  need  to  align  the  two  optical  systems  to 
each  other  with  subarcsecond  accuracy.  More  important,  the 
relative  alignment  of  two  separate  optical  systems  would  drift  with 
time  and  temperature,  whereas  the  star  trackers  for  the  Mark  III 
interferometer  will  maintain  alignment  indefinitely. 

The  photon  camera  used  in  the  star  tracker  is  a  commercial 
detector  based  on  the  PAPA  camera  (Papaliolios  et  al.,  1985)  with 
256  x  256  pixels.  The  detector  produces  a  pulse  with  an  x,  y 
address  for  every  detected  photon.  The  digital  output  of  the 
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Fig.  2.  Optical  schematic  of  the  Mark  III  interferometer 


camera  is  interfaced  to  a  high  speed  signal  processor  which  locates 
the  centroids  of  the  two  stellar  images.  Two  additional  optical/ 
mechanical  components  are  shown  in  Fig.  3.  The  variable  neutral 
density  filter  is  used  to  reduce  the  photon  flux  to  the  camera  to 
prevent  overloading  of  the  camera  electronics  when  observing 
bright  stars.  The  setting  of  the  filter  is  controlled  by  the  star  tracker 
computer.  There  are  also  two  shutters,  one  in  each  arm  of  the 
interferometer,  which  are  used  for  star  acquisition.  When  the 
siderostats  stop  slewing,  the  two  stellar  images  may  be  sufficiently 
far  from  their  correct  positions  that  proper  identification,  i.e., 
which  star  came  from  which  siderostat,  may  not  be  possible.  The 
shutters  are  used  to  remove  this  ambiguity. 

The  camera  is  the  error  sensor.  The  actuator  is  a  PZT-driven 
gymbaled  mirror.  This  5-inch  mirror  is  driven  in  two  orthogonal 
directions  with  a  30"  peak-to-peak  motion.  The  bearings  are  cross 
flexure  pivots  and  the  mirror  mount  has  its  lowest  resonance  at 
200  Hz. 

The  heart  of  the  star  tracker  is  the  control  system;  it  makes 
extensive  use  of  optimal  estimation  techniques  and  only  a  brief 
outline  is  given  here.  A  detailed  description  is  given  by  Clark  et  al. 
(1986).  The  control  problem  is  separated  into  a  spatial  filtering 
part  and  a  temporal  filtering  part.  The  spatial  filter  takes  the 
photon  addresses  generated  by  the  photon  camera  and  derives  the 


centroid  of  the  star.  The  temporal  filter  takes  current  and  past 
centroid  estimates  and  generates  a  set  of  voltages  to  drive  the 
piezoelectric  actuators. 

The  image  of  the  star  on  the  camera  face  is  the  diffraction 
pattern  of  an  annular  telescope  aperture.  A  cross-correlation 
centroid  estimator  was  implemented.  The  estimated  x  coordinate 
of  the  centroid  is  given  by 

x=$f(x,y)hx(x,y)dxdy,  (1) 

where /  is  the  detected  image  and  hx  is  the  correlation  function  for 
the  x  direction.  The  noise  in  the  detected  image  is  photon  noise 
with  a  variance  given  by 

o2f=f(x,y)-  (2) 


Since  the  theoretical  shape  of  the  image  is  known,  there  is  a  closed 
form  solution  for  the  optimal  correlation  function,  i.e.,  that 
function  h(x,y)  which  gives  a  minimum  variance  for  x  due  to 
photon  noise  in  the  measurement  of  f(x,y).  The  optimal 
correlation  function  is  given  by 


hopt  (“>$  = 


[~df(<*,P)ldx] 

f(«J) 


(3) 
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OPTICAL  LAYOUT  OF  STAR  TRACKER 


Fig.  3.  Optical  schematic  of  star  tracker 


Table  1 


Correlation  function 

Single  photon  variance 

(arcsec)2 

Optimal  [Eq.  (3)] 

0.303 

Quad  cell  [Eq.  (4)] 

3.09 

Linear  [Eq.  (5)] 

6.14 

This  type  of  centroid  estimator  has  limited  dynamic  range  but  very 
low  variance.  An  analysis  of  this  estimator  for  our  annular 
aperture  image  was  made  relative  to  centroid  estimators  using 
quad  cell  and  linear  correlation  functions: 

Quad  cell:  h=_y  x<Q,  (4) 

linear:  h  =  x.  (5) 

The  results  are  shown  in  Table  1.  Note  that  the  variance  for  the 
optimal  correlator  estimator  is  about  an  order  of  magnitude  better 
than  a  quad  cell. 

The  correlation  integral  [Eq.  (1)]  is  performed  by  a  custom 
digital  signal  processor.  The  processor  calculates  the  integral  on  a 
photon  by  photon  basis,  rather  than  a  pixel  by  pixel  basis.  This  is 
more  efficient  computationally  when  the  number  of  detected 
photons  is  much  less  than  the  number  of  pixels  in  the  image. 
Currently,  centroid  estimates  are  calculated  every  10  ms.  The 
signal  processor  can  operate  at  photon  rates  well  above  the 
capabilities  of  the  camera.  The  optical  correlation  functions,  one 
in  x,  and  one  in  y  for  each  of  the  two  trackers  are  downloaded  to 
the  correlation  processor  from  the  main  control  computer.  The 


calculation  of  the  optimal  correlation  function  is  not  done  in  real 
time. 

The  second  part  of  the  control  system  is  the  temporal  filter. 
While  the  spatial  filter  is  implemented  mostly  in  digital  hardware 
and  microcode,  the  temporal  filter  is  implemented  entirely  in 
software  on  the  star  tracker  control  computer.  The  temporal  filter 
is  a  discrete  time  Kalman  filter.  The  Kalman  filter  assumes  a  state 
variable  formulation  of  the  tracking  problem  and  is  the  linear 
filter  which  minimizes  the  mean  square  tracking  error.  This  filter 
takes  into  account  the  observation  noise  (photon  noise)  by 
adaptively  varying  the  tracker  bandwidth  according  to  the  noise 
variance.  The  output  of  the  filter  is  scaled  and  then  directed  to  a 
D/A  converter  and  high  voltage  amplifier  to  drive  the  PZTs. 

In  addition  to  the  software  for  the  Kalman  filter,  there  is 
software  for  acquisition,  data  recording,  and  for  utilities  such  as 
moving  the  shutters  and  controlling  the  camera.  There  is  also 
software  for  complete  automatic  operation  and  for  off-line 
diagnostics  and  analysis  of  tracker  performance. 

2.3.  Optical  delay  line 

Two  optical  delay  lines  are  used  to  equalize  the  paths  from  the  star 
to  the  beamsplitter  via  the  two  siderostats  12  m  apart.  The 
mechanical  and  optical  arrangement  is  shown  in  Fig.  4.  The 
mechanical  design  is  similar  to  that  used  in  delay  lines  for  high 
resolution  Fourier  spectroscopy  (Connes,  1975).  The  optical  delay 
is  changed  by  moving  cat’s  eye  retroreflectors  along  tracks  which 
are  housed  in  two  25-ft  vacuum  chambers.  The  peak-to-peak 
optical  path  difference  of  this  system  is  20  m,  and  the  delay  lines 
can  be  positioned  under  laser  control  to  an  accuracy  of  10  nm,  for 
a  dynamic  range  of  2  109  resolution  elements  in  optical  path 
difference.  The  maximum  slew  rate  is  1  m  per  second  of  optical 
path  difference.  Twenty  seconds  is  required  to  go  from  one  end  to 
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Fig.  4.  Three  stage  (PZT,  Voice  Coil,  Motor)  optical  delay  line 


Fig.  5.  Heterodyne  laser  interferometer  for  delay  line 


the  other.  The  maximum  tracking  rate  is  2  mm  s  “  \  at  which  point 
the  rms  position  accuracy  has  degraded  to  20  nm. 

Figure  4  shows  the  optical  path  of  the  laser  beam  that 
measures  the  position  of  the  cat’s  eye  retroreflector.  The  cat’s  eye 
is  an  8-inch  f/3  parabola  with  a  1-inch  flat  mirror  at  its  focus.  The 
1-inch  mirror  is  mounted  on  a  PZT  stack  that  has  a  5-pm  peak-to- 
peak  range.  The  cat’s  eye  is  suspended  on  a  flexure  bearing  over  a 
belt  driven  translation  stage.  A  voice  coil  positions  the  cat’s  eye 
over  the  stage,  while  a  motor  at  one  end  of  the  25-ft  track  moves 
the  stage  via  a  timing  belt.  The  voice  coil  has  a  2-mm  range  and  the 
motor  has  a  9-pm  step  size.  The  three  actuators  -  PZT,  voice  coil, 
and  motor  -  have  considerable  overlap  in  terms  of  resolution  and 
dynamic  range. 

The  delay  line  sensor  is  a  heterodyne  laser  interferometer 
(Fig.  5).  The  light  from  a  stabilized  laser  is  split  into  two  beams  by 
a  polarizing  beamsplitter.  They  are  then  frequency  shifted  by  two 
acousto-optic  bragg  cells  with  a  2-MHz  frequency  difference. 
When  recombined,  the  beam  has  its  two  orthogonal  polarizations 
at  different  frequencies  with  a  2-MHz  frequency  difference.  The 
laser  is  a  commercial  transverse  Zeeman  stabilized  HeNe  laser 
with  a  one-day  stability  of  200  kHz  at  line  center,  474  THz.  The 


interferometer  is  a  standard  polarization  heterodyne  interfero¬ 
meter.  The  phase  of  the  2-MHz  signal  from  the  interferometer  is 
compared  to  the  2-MHz  reference  (Fig.  5)  and  a  one  wavelength 
change  in  the  monitored  optical  path  causes  a  2  n  change  in  the 
relative  phases  of  the  interferometer  and  reference  2-MHz  signals. 
A  custom  digital  phase  comparator  circuit  counts  the  fringes  with 
2/64  resolution  at  a  maximum  slew  rate  of  2  10 6  fringes  per 
second. 

Since  the  desired  optical  path  difference  is  given  by  the 
difference  between  the  positions  of  the  two  delay  lines,  it  is  only 
necessary  to  actively  control  one  delay  line.  Thus,  vibrations  in 
one  delay  line  are  compensated  by  driving  the  other  in  exactly  the 
same  way. 

The  laser  interferometers  are  read  at  a  1-kHz  rate  by  the  fringe 
tracker  computer.  The  difference  between  the  two  interferometers 
is  the  measured  position,  which  is  subtracted  from  the  desired 
position  to  give  an  error  signal.  The  PZT  servo  is  a  simple  first- 
order  system,  i.e.,  the  PZT  drive  voltage  is  just  the  integral  of  the 
error  signal.  The  necessary  integral  is  performed  in  digital 
hardware,  and  the  velocity  of  the  PZT  mirror  is  updated  every 
millisecond  by  the  computer.  The  closed-loop  bandwidth  of  this 
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1  servo  is  approximately  100  Hz,  for  a  time  constant  of  1  to  2  ms. 
i  The  small-signal  bandwidth  of  the  entire  delay  line  is  the 
!  bandwidth  of  this  innermost  servo. 

[  The  voice  coil  servo  is  designed  to  keep  the  PZT  system  within 
[  its  5-pm  dynamic  range.  The  cat’s  eye  flexure  suspension  has  a  1 . 5- 
i  Hz  resonance.  Since  a  5-pm  accuracy  requires  a  response  signifi- 
i  cantly  faster  than  1.5  Hz,  the  voice  coil  servo  has  to  be  at  least  a 
second-order  servo.  The  voice  coil  servo  is  implemented  entirely  in 
software.  The  error  signal  input  to  the  voice  coil  servo  is  the 
distance  the  PZT  mirror  is  from  the  center  of  its  range.  The  output 
is  a  current  to  the  voice  coil.  The  closed-loop  bandwidth  of  this 
servo  is  approximately  10  Hz. 

The  motor  servo  is  designed  to  keep  the  voice  coil  in  the  center 
of  its  range.  The  input  to  this  servo  is  the  position  of  the  voice  coil 
as  measured  with  a  lateral-cell  transducer.  The  motor  is  a 
microstep-driven  stepper  motor,  interfaced  to  the  computer 
through  a  programmable  pulse  generator.  In  the  tracking  mode, 
it  operates  as  a  first-order  servo,  with  a  small-signal  closed-loop 
bandwidth  of  about  1  Hz. 

The  delay  line  control  system  was  designed  to  interact  with  the 
rest  of  the  interferometer  through  a  simple  interface.  The  rest  of 
the  interferometer  does  not  see  three  separate  actuators  and 
control  loops  and  a  laser  interferometer.  Rather,  the  delay  line 
subsystem  is  programmed  by  setting  its  position  and  velocity.  At  a 
given  time,  the  delay  line  is  programmed  to  move  to  a  given 
position.  Thereafter,  the  optical  delay  is  changed  at  a  rate 
determined  by  the  programmed  velocity.  The  delay  line  will  move 
at  this  rate  until  the  programmed  position  or  velocity  is  changed. 

The  composite  three-stage  servo  system  has  a  moderately  high 
bandwidth  when  the  error  signal  is  small  (i.e.,  when  delay  line 
position  minus  target  position  is  less  than  about  4  pm).  However, 
when  switching  from  one  star  to  the  next,  the  programmed 
position  can  suddenly  change  by  15  m.  Since  the  PZT  and  voice 
coil  have  small  dynamic  ranges,  and  in  addition  have  limits  on 
their  velocity  and  acceleration,  these  servo  systems  will  not  behave 
properly  in  response  to  such  a  large  step.  Left  uncorrected,  the 
resultant  oscillations  are  sufficiently  violent  to  misalign  optics. 
Therefore,  in  addition  to  the  linear  tracking  mode  of  operation 
described  above,  the  delay  line  has  several  other  modes  of 
operation.  A  significant  portion  of  the  delay  line  servo  code  detects 
input  conditions  which  would  require  a  change  in  the  servo  mode 
so  that  the  delay  line  remains  well  behaved  under  all  conditions. 

Table  2  shows  the  various  modes  of  operation  of  the  three 
delay-line  loops.  The  modes  of  the  various  servos  depend  on  the 
inputs  to  the  delay  line.  Rapid  changes  in  the  programmed 
position  or  velocity  is  the  usual  reason  for  changes  in  operating 
mode.  Hysteresis  in  the  transition  thresholds  prevents  oscillations 
caused  by  changes  in  the  servo  parameters  between  modes.  Inputs 
from  the  limit  switches  can  also  cause  changes  in  the  servo  modes. 
For  example,  two  limit  switches  are  used  at  each  end  of  the  track. 
The  first  informs  the  computer  that  the  delay  line  should  start 
decelerating  in  order  to  prevent  the  delay  line  from  moving 
beyond  the  second  limit  switch,  which  disables  the  motors.  Other 
interlocks  prevent  the  laser  fringe  counter  from  running  at  2 
million  fringes  per  second,  taking  the  delay  line  with  it,  if  the  laser 
becomes  misaligned.  For  a  more  detailed  description  of  the  delay 
line,  the  reader  is  referred  to  Holm  (1986). 

2.4.  Stellar  fringe  detector 

With  proper  operation  of  the  star  tracker  and  delay  line,  the  two 
beams  combined  at  the  beamsplitter  should  interfere.  This  section 


Table  2.  Inputs:  target  position,  target  velocity,  current  laser 
positions  (2),  limit  switches  (8),  voice  coil  position 


Servo 

Modes 

PZT 

Linear  track 

Off 

Voice  coil 

Linear  track 

Damp 

Off 

Motor 

Linear  track 

Accel./decel. 

Constant  vel.  Off 

Fig.  6.  Stellar  fringe  modulation/demodulation 


describes  the  procedure  by  which  the  amplitude  and  phase  of  the 
interference  pattern  are  measured. 

In  the  Mark  III  interferometer,  the  interfering  wavefronts  are 
parallel.  Hence,  the  two  outputs  of  the  beamsplitter  are  uniform 
over  the  aperture  and  complementary.  If  the  phase  difference  is  90 
degrees,  all  then  light  will  go  to  one  of  the  two  outputs;  if  the  phase 
difference  is  —90  degrees,  the  light  will  go  to  the  other  output. 

The  fringe  detection  algorithm  has  been  described  previously 
(Shao  and  Staelin,  1977).  The  basic  technique  uses  path  length 
modulation  and  synchronous  demodulation.  The  PZT  controlled 
mirror  in  one  of  the  vacuum  delay  lines  is  modulated  with  a  500- 
Hz  triangle  wave.  If  the  path  difference  in  the  two  arms  of  the 
interferometer  is  within  the  coherence  length  of  the  light  (a  few 
wavelengths),  the  intensity  of  light  at  the  output  of  the  beamsplit¬ 
ter  will  be  as  shown  in  Fig.  6.  During  one  sweep  of  the  PZT,  the 
detected  photons  are  grouped  into  four  time  bins  denoted  by  A,  B, 
C,  and  D.  The  amplitude  and  phase  of  the  fringe  are  given  by 


4>  =  tanI(, C-A )  4’ 

(6) 

Amp  =  j/(C-^)2  +  (Z>-5)2 . 

(7) 

In  the  above  equations,  the  C— A  and  D  —  B  terms  are  the  sine  and 
cosine  components  of  the  fringe. 

The  measured  phase  is  unwrapped  so  that  if  the  atmosphere 
moves  the  fringe  position  beyond  +  n,  the  correct  error  signal  will 
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!  be  used  to  track  the  fringe.  The  unwrapped  phase  is  used  as  the 
!  error  signal  in  a  first-order  servo.  As  mentioned  before,  the  delay 
[  line  is  programmed  by  position  and  velocity.  The  unwrapped 
[  phase  times  a  gain  term  is  added  to  the  previous  programmed 
>  position  of  the  delay  line  to  yield  the  new  programmed  position. 

!  The  programmed  velocity  is  the  sidereal  rate.  The  closed-loop 
bandwidth  of  this  servo  is  about  20  Hz.  A  similar  fringe  tracker 
was  used  on  the  Mark  I  and  Mark  II  interferometers  to  track 
stellar  fringes. 

The  fringe  tracking  servo,  implemented  entirely  in  software,  is 
the  simplest  part  of  the  stellar  fringe  detection  subsystem.  Other 
software  components  of  the  fringe  detection  subsystem  perform 
functions  such  as  stellar  fringe  search  and  initial  acquisition, 
fringe  recovery  after  tracker  failure,  central  fringe  identification, 
and  automatic  transition  between  operating  modes.  In  addition  to 
the  above  software  components,  other  system  components  that 
affect  the  operation  of  the  stellar  fringe  detector  are  the  PZT  stroke 
servo,  and  the  optics  for  two  color  fringe  measurements  and  for 
accurate  visibility  measurements. 

The  instrument  was  designed  to  make  accurate  astrometric 
measurements  and  accurate  stellar  diameter  measurements  on  as 
many  stars  as  possible.  The  design  of  a  stellar  fringe  tracker  is 
complicated  by  the  fact  that  the  requirements  of  high  tracker 
sensitivity,  low  systematic  astrometric  error,  and  low  visibility 
errors  are  not  entirely  consistent. 

The  stroke  of  the  PZT  can  be  varied  by  one  of  the  control 
computers.  Since  the  piezoelectric  constant  is  temperature  de¬ 
pendent,  the  peak-to-peak  amplitude  of  the  PZT  triangle  wave  is 
monitored  with  a  laser  interferometer  and  the  drive  voltage  is 
adjusted  to  maintain  the  proper  amplitude  to  within  0.5  nm.  For 
stellar  diameter  measurements,  the  stroke  is  set  at  the  effective 
wavelength  of  the  narrow  bandpass  filter.  For  astrometry,  the 
stroke  is  set  at  the  effective  wavelength  of  the  red  channel. 

When  the  interferometer  is  set  to  point  at  a  star,  the  delay  line 
is  programmed  to  go  to  the  expected  position  of  the  fringe  plus 
60  pm.  From  past  experience,  the  calculated  fringe  position  is 
never  more  than  30  pm  in  error  if  the  baseline  vector  used  in  the 
calculation  is  derived  from  stellar  fringe  data  less  than  a  week  old. 
The  delay  line  is  programmed  to  scan  for  the  fringe  at  approxi¬ 
mately  60  pm  per  second.  The  fringe  amplitude  is  calculated  every 
atmospheric  coherence  time,  and  if  a  fringe  is  detected,  the  scan  is 
stopped  and  the  servo  turned  on. 

When  the  tracker  loses  the  fringe,  the  delay  line  is  moved  back 
5  pm  and  a  search  is  restarted.  After  initial  fringe  acquisition,  the 
search  space  is  restricted  to  a  20-pm  interval  centered  about  the 
most  recent  position  where  the  fringe  was  found.  All  of  these 
search  parameters  are  superimposed  on  the  fringe  motion  due  to 
earth  rotation,  which  is  calculated  using  a  baseline  vector  typically 
derived  from  the  previous  night’s  data.  When  the  seeing  is 
exceptionally  good,  i.e.,  sub-arc-second  seeing,  the  fringe  tracker 
will  stay  locked  on  one  fringe  for  many  minutes,  for  up  to  105 
atmospheric  coherence  times. 

If  the  fringe  stays  locked  for  more  than  one  second,  a  central 
fringe  identification  algorithm  is  enabled.  The  visibility  during 
that  one-second  interval  is  compared  with  a  threshold.  If  the 
visibility  is  below  the  threshold,  it  is  assumed  that  the  tracker  is  on 
the  wrong  fringe  and  the  tracker  will  jump  one  fringe  to  the  right 
or  left.  The  direction  of  the  tracker  jumps  is  determined  by  the 
visibility  gradient.  Since  the  tracker  has  an  rms  phase  error  of 
about  1/6  of  a  fringe,  the  visibility  gradient  can  be  determined  by 
looking  at  the  distribution  of  visibility  versus  phase.  Since  the 
central  fringe  has  the  highest  visibility,  the  tracker  will  usually  hop 
in  that  direction. 


The  operation  of  the  various  systems  which  comprise  the  star 
light  fringe  detector  requires  that  a  number  of  quantities:  phase, 
amplitude,  SNR,  visibility  gradient,  etc.  be  calculated  in  real  time. 
Although  the  system  is  rather  complicated,  the  resulting  tracker  is 
extremely  robust.  The  interferometer  can  be  commanded  to  point 
at  a  star  and  the  tracker  will  lock  onto  the  fringe  within  90  s.  Thus, 
the  software  complexity  was  introduced  to  enable  the  instrument 
to  operate  efficiently  under  a  wide  variety  of  conditions  with  a 
minimum  of  operator  intervention. 

2.5.  Laser  metrology 

During  initial  operation  in  1986,  the  siderostat  laser  metrology 
system  was  not  operational.  However,  the  laser  system  is  impor¬ 
tant  for  astrometry  and  should  become  operational  in  the  summer 
of  1987.  The  basic  laser  interferometer  used  to  measure  changes  in 
instrument  geometry  was  described  in  the  delay  line  section  of  this 
paper.  Two  such  interferometers  measure  the  position  of  the  cat’s 
eye  retro  reflectors.  An  additional  12  interferometers  will  be  used 
to  measure  the  position  of  the  three  siderostat  mirrors. 

The  metrology  system  measures  the  motion  of  the  pivot  point 
of  the  siderostat  by  measuring  the  distance  between  four  reference 
corner  cubes  in  an  invar  plate,  which  is  bonded  to  the  concrete 
pedestal,  and  the  surface  of  a  zerodur  sphere  attached  to  the  mirror 
as  shown  in  Fig.  7.  The  sphere  is  made  accurate  to  2/10  and  the 
center  of  the  sphere  is  coincident  with  the  mechanical  pivot  point 
to  20  pm.  The  interferometers  have  5  nm  resolution;  the  desired 
accuracy  of  the  pivot  point  monitor  is  100  nm  for  one-milliarcsec 
astrometry  with  a  20-m  baseline.  The  measurements  from  these 
lasers  will  be  used  to  correct  for  mechanical  errors  in  the  siderostat 
as  well  as  thermal  motion  of  the  siderostat,  as  described  below. 

The  basic  model  for  the  instrument  is  given  by 

Fringe  Pos  =  S  •  B  -f  C ,  (8) 

S  is  the  unit  vector  to  the  star, 

B  is  the  baseline  vector, 

C  is  the  delay  offset. 

If  the  siderostats  were  perfect,  the  baseline  vector  would  be  the 
vector  that  joined  the  pivot  points  of  the  two  siderostats.  If  the 
delay  lines  were  commanded  to  move  to  their  zero  positions,  the 
delay  offset  C  would  be  the  difference  in  the  optical  paths  from  the 
beamsplitter  to  the  pivot  points  of  the  two  siderostats.  Motion  of 
the  siderostat  pivot  points  affects  both  B  and  C.  Motion  of  the 
internal  optics  only  affects  C. 

With  our  12-m  baseline  on  Mt.  Wilson,  the  rms  atmospheric 
fringe  motion  was  6  pm  under  conditions  of  good  seeing.  For 
stellar  diameter  measurements  with  the  fringe  tracker  running, 
neither  instrument  drifts  nor  atmospheric  motion  is  important. 
However,  for  astrometric  measurements,  the  measured  fringe 
position  in  Eq.  (8)  is  made  relative  to  the  baseline  vector  and  delay 
offset.  Thus,  it  is  necessary  that  both  B  and  C  be  constant  over  a 
period  of  time  that  is  sufficiently  long  for  the  interferometer  to  be 
able  to  observe  the  fringe  position  of  a  number  of  stars. 

The  parameters  B  and  C  can  change  for  two  reasons:  imperfect 
siderostat  bearings  and  thermal  drift.  Ideally,  the  siderostat  pivots 
about  a  point  on  the  surface  of  the  mirror.  In  practice,  the  mirror 
translates  as  it  rotates  in  a  rather  unpredictable  way  due  to 
irregularities  in  the  bearings  and  flexure  of  the  mirror  cell.  The 
laser  metrology  system  is  designed  to  measure  these  imperfections. 
The  result,  in  theory,  is  that  baseline  drifts  will  be  limited  to 
thermal  drift  of  the  top  of  the  concrete  pier  to  which  the  reference 
corner  cubes  in  the  laser  system  are  attached.  With  rapid  switching 
between  stars,  only  baseline  drift  with  a  time  scale  faster  than  the 
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Fig.  7.  Laser  metrology  system  for  baseline  vector  monitor 


SIDEROSTAT  WITH  4-LASER  PLATE 


switching  cycle  results  in  astrometric  errors  for  large  angle  relative 
astrometry. 

In  order  to  monitor  the  position  of  the  pivot  point,  a  minimum 
of  three  interferometers  is  necessary.  The  fourth  interferometer 
provides  the  redundancy  necessary  to  solve  for  the  position  of  the 
reference  corner  cubes.  The  interferometers  measure  the  position 
of  the  center  of  the  sphere.  Certain  parameters,  such  as  the 
distance  between  the  center  of  the  sphere  and  the  surface  of  the 
mirror  must  be  derived  from  stellar  data. 

The  mechanical  accuracy  of  the  instrument  without  the  laser 
metrology  system  is  approximately  10-15  pm  if  the  baseline  was 


calculated  from  the  previous  night’s  stellar  data.  The  size  of  this 
error  implies  that  searching  for  the  fringe  would  take  only  a  few 
hundred  milliseconds,  an  insignificant  amount  of  time  compared 
to  the  30  to  60  s  needed  for  the  siderostats  to  slew  to  a  new  star. 
The  mechanical  errors  are  roughly  equivalent  to  the  errors  due  to 
atmospheric  turbulence,  or  to  the  errors  in  the  primary  mirror  of  a 
high  quality  large  astronomical  telescope.  With  the  laser  system 
we  hope  to  improve  this  accuracy  by  two  orders  of  magnitude. 

While  this  higher  mechanical  accuracy  is  not  necessary  for 
stellar  diameter  measurements,  it  is  crucial  for  ultra-high  accuracy 
astrometry.  This  level  of  accuracy  is  also  necessary  for  future 
space-based  phase  coherent  long  baseline  interferometers. 


2.6.  Siderostat  control 

The  siderostat  control  subsystem  is  in  many  ways  similar  to  the 
computer  controlled  pointing  system  for  a  modern  astronomical 
telescope.  The  siderostat  itself  is  a  two  axis  motorized  gymbaled 
mirror  mount  (Fig.  8),  with  the  fixed  axis  tilted  75°  from  vertical 
towards  the  West.  The  siderostats  use  conventional  precision  worm 
and  wheel  gears  and  microstepped  stepper  motors.  The  control 
computer  sets  programmable  pulse  generators  to  control  motor 
velocity.  The  average  step  size  is  0"02  to  provide  smooth  tracking, 
and  the  maximum  slew  rate  is  1?8  per  second.  Typical  backlash  in 
the  gear  train  is  less  than  1". 

The  control  system  performs  several  tasks.  One  is  open  loop  or 
blind  pointing,  where  a  geometric  model  of  the  siderostat  is  used 
to  point  the  siderostat  and  track  the  star.  A  second  task  is  closed 
loop  pointing,  where  the  siderostat  is  moved  to  keep  the  high¬ 
speed  PZT  tilt  mirrors  in  the  center  of  their  range.  The  third  task  is 
updating  the  parameters  of  the  geometric  model  when  a  star  is 
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!  Table  3.  Siderostat  model  parameters 


[ 

I  Az,  E1(F) 
_ f) 

u  axis 
fj 

^mirror 


Angular  orientation  of  feed  beam  (internal 
alignment) 

Azimuth  and  elevation  of  fixed  axis 
Angle  between  fixed  and  movable  axis 
Angle  between  movable  axis  andmirror  normal 
Zero  point  of  motor  position  on  the  two  axes 


acquired.  In  addition,  the  siderostat  control  system  controls  the 
sequencing  of  observations  during  the  night. 

The  geometry  of  the  siderostat  described  earlier  was  that  of  an 
ideal  siderostat.  However,  the  stationary  axis  may  not  be  exactly 
1 5  degrees  from  horizontal;  it  may  not  be  pointed  due  west,  and  the 
two  axes  may  not  be  exactly  perpendicular.  Altogether,  there  are 
eight  parameters  in  the  geometric  model  of  the  siderostat.  With 
values  of  these  eight  parameters,  there  is  a  closed-form  relation 
between  the  angular  positions  of  the  motors  and  the  pointing  of 
the  instrument.  It  is  also  possible  to  invert  this  problem  in  closed 
form.  By  pointing  the  siderostat  at  a  number  ( >  4)  of  bright  stars 
whose  coordinates  are  known,  and  recording  the  positions  of  the 
motors  when  the  stars  are  centered  in  the  field  of  view  of  the 
interferometer,  the  eight  parameters  of  the  model  can  be  deter¬ 
mined  in  a  least  squares  fashion.  In  our  initial  setup,  we  used  the 
FK4  coordinates  of  a  half  dozen  stars  precessed  to  the  night  of 
observation. 

Since  the  siderostats  direct  the  starlight  into  the  central  beam¬ 
combining  building,  the  alignment  of  the  interferometer  optics 
inside  the  building  will  affect  the  pointing  of  the  siderostat.  The 
internal  alignment  of  the  interferometer  is  represented  by  two  of 
the  eight  parameters  describing  the  siderostat.  All  of  the  para¬ 
meters  are  listed  in  Table  3.  No  attempt  was  made  to  model  errors 
in  the  precision  worm  and  wheel  gear  or  to  model  flexure.  It  is 
most  likely  that  the  residual  errors  are  primarily  in  the  gears. 
Residuals  from  a  least-squares  fit  for  the  siderostat  parameters  are 
of  order  5".  Typical  pointing  errors  for  stars  which  are  not  used  in 
the  least-squares  fit  are  approximately  7"  to  10". 

However,  each  new  star  that  is  observed  provides  some 
information  on  the  eight  siderostat  parameters.  Hence,  rather 
than  an  off-line  least-squares  solution,  a  real  time  algorithm  could 
be  implemented  to  update  the  eight  parameters  every  time  the 
instrument  locked  onto  a  new  star.  In  the  Mark  III,  a  Kalman 
filter  is  used  to  update  the  parameter  vector.  In  operation,  the 
initial  values  of  the  parameters  at  the  beginning  of  the  night  are  the 
parameters  derived  from  the  previous  night’s  observations.  These 
parameters  are  updated  for  every  new  star  and,  under  normal 
conditions,  the  initial  values  would  be  essentially  lost  after 
approximately  6  stars.  Thereafter,  the  parameter  estimates  would 
reflect  thermal  drift  of  the  instrument  during  the  night. 

This  complex  siderostat  model  is  used  only  for  initial  acqui¬ 
sitions.  After  the  star  tracker  locks  onto  the  star,  the  positions  of 
the  PZT  tilt  mirrors  are  sent  to  the  siderostat  control  computer. 
The  siderostat  drive  rates  for  the  motors  are  then  modified  to  keep 
the  PZTs  in  the  center  of  their  range.  Since  the  PZT  tilt  mirrors  in 
effect  represent  realignment  of  the  internal  optics,  the  corrections 
to  the  sideral  rate  are  generated  by  changing  the  parameters 
^F  5  4^F  * 

The  siderostat  control  computer  also  calculates  the  theoretical 
fringe  position  and  velocity  for  the  optical  delay  line.  These  values 


are  sent  to  the  delay  line  control  computer  through  a  shared 
multiport  memory  every  2  s.  The  theoretical  fringe  position  is 
calculated  from  a  baseline  vector  derived  from  the  previous  night’s 
data,  in  the  future,  we  plan  to  update  the  baseline  vector  with  a 
Kalman  filter  in  the  same  way  the  siderostat  parameters  are 
updated.  On  the  first  night  of  interferometric  observations,  the 
baseline  vector  was  determined  by  conventional  surveying  techni¬ 
ques.  The  surveying  error  was  approximately  1  mm,  or  about  20" 
of  baseline  tilt. 

The  last  function  of  the  siderostat  control  computer  is  the 
sequencing  of  observations.  An  abbreviated  version  of  the  FK4 
catalog  is  kept  on-line  (In  1987,  the  FK5  catalog  will  be  kept  on¬ 
line).  A  menu  driven  program  enables  the  operator  to  list  and  plot 
the  stars  visible  at  the  present  time  and  to  easily  define  the  stellar 
observation  sequence.  The  siderostat  control  computer  also 
communicates  with  the  other  two  control  computers  to  enable 
fully  automatic  operation  once  the  observing  sequence  is  pro¬ 
grammed.  The  siderostat  control  computer  signals  the  star  tracker 
computer  to  initiate  a  star  acquisition  sequence.  When  the  star 
tracker  has  locked  onto  a  star  it  signals  the  fringe  tracker 
computer  to  initiate  a  fringe  acquisition  sequence.  When  the 
allotted  observing  time  for  a  star  is  over,  the  siderostat  control 
computer  signals  the  other  two  computers  to  disengage  their  servo 
systems.  For  additional  information,  the  reader  is  referred  to 
Kelly  (1986). 


3.  Operational  characteristics 

The  instrument  was  built  and  debugged  over  a  33  month  period. 
Fringes  were  observed  on  the  first  night  of  operation,  September 
12,  1986.  Over  the  following  months,  observational  techniques 
were  refined  and  automated.  Figure  9  is  a  plot  of  the  observations 
made  on  1986,  November  11.  The  vertical  axis  is  the  delay  line 
position  where  fringes  were  observed.  Each  dot  represents  a  100-s 
observation.  A  total  of  165  observations  of  23  stars  were  made 
over  a  10.5-h  period.  Stars  over  a  declination  range  of  68°  and  an 
hour  angle  range  of  8  h  were  observed.  All  these  observations  were 
of  stars  3.0  magnitude  or  brighter. 

The  initial  measurements  in  1986  will  be  discussed  in  a  series  of 
separate  papers.  Wide-angle  one-color  astrometry  shows  night-to¬ 
night  repeatability  of  about  0"04  (Mozurkewich  et  al.,  1987),  and 
the  formal  error  for  seven  nights  of  data  is  0"02  in  declination. 
Significant  improvements  are  anticipated  when  two-color  analysis 
is  applied  and  when  the  laser  siderostat  monitor  is  installed  in  the 
summer  of  1987.  Preliminary  measurements  on  the  atmospheric 
limitations  of  two  color  absolute  astrometric  measurements 
indicate  that  single  observations  of  ~  100  s  should  yield  accuracies 
of  about  0"01,  a  factor  of  5  improvement  over  one-color 
observations  (Colavity  et  al.,  1987).  Two-color  observations 
should  approach  0"003  for  a  single  night’s  observation.  Last  of  all, 
we  have  made  preliminary  stellar  diameter  measurements  of 
/?Gem  and  aCMi.  Assuming  a  uniform  disk,  we  measure 
6.15  +  0.25  mas  for  aCMi  and  8.70  +  0.15  mas  for  /?Gem.  The 
uniform  disk  diameter  for  aCMi  is  not  in  agreement  with  prior 
measurements  (Hanbury  Brown  et  al.,  1974)  in  part  due  to  the 
difference  in  observation  wavelength  and  uniform  disk  assump¬ 
tion.  A  much  more  careful  discussion  of  accurate  stellar  diameter 
measurements  is  given  elsewhere  (Shao  et  al.,  1988). 

The  instrument  was  never  operated  near  the  limiting  magni¬ 
tude  of  the  fringe  tracker  because  the  temporary  star  tracker  we 
were  using  could  not  track  stars  fainter  than  3  mag.  The  tracker 
using  the  photon  camera  should  be  capable  of  operation  on  faint 
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Fig.  9.  Position  of  delay  line  where  fringes  were  observed,  165  observations  of  23  stars 


sources,  to  12  mag  under  good  seeing  conditions  when  only 
mechanical  errors  need  to  be  tracked.  When  the  new  star  tracker  is 
installed  in  the  spring  of  1987,  the  limiting  magnitude  will  be  set  by 
the  fringe  detection  system.  We  plan  to  operate  the  interferometer 
in  two  distinct  modes  that  have  two  very  different  magnitude 
limits.  The  two  modes  represent  two  types  of  interferometers.  One 
is  a  closed  loop  phased  interferometer.  The  second  is  an  open  loop 
coherent  interferometer  or  absolute  interferometer.  We  plan  to 
test  this  second  mode  of  operation  in  the  next  two  years. 

3.1.  Phase-coherent  interferometry 

The  Mark  III  was  designed  to  be  an  active  phased  interferometer. 
In  this  mode  of  operation,  the  fringe  tracker  compensates  for  all 
optical  path  fluctuations,  including  atmospheric  turbulence  to 
~/l/6rms.  Fringe  phase  measurements  in  two  colors  are  used  for 
astrometric  purposes,  while  light  from  a  bandpass  filter,  <  400  A, 
is  used  for  stellar  diameter  measurements.  It  is  also  possible  to 
measure  the  diameter  of  a  star  at  emission  or  absorption  lines.  For 
diameter  measurements  in  narrow  spectral  lines,  we  make  use  of 
the  phase  coherent  property  of  an  active  interferometer.  For 
narrow  spectral  lines,  the  photon  flux  in  the  line  might  be  10“ 3  of 
the  total  flux  and  the  number  of  photons  per  atmospheric 
coherence  time  may  be  as  low  as  10“ 2 .  However,  since  the  active 
optics  system  is  phase  coherent,  the  coherent  integration  time  is 
not  limited  to  the  atmospheric  coherence  time.  By  choosing  a 
longer  coherent  integration  time,  >  1  s,  a  fringe  visibility  measure¬ 
ment  can  be  made  in  104s  with  a  signal-to-noise  ratio  of  100. 

For  astrometry,  the  fringe  phase  and  the  position  of  the  optical 
delay  line  is  recorded  when  the  fringe  tracker  is  locked  on  the 


central  fringe.  In  this  case,  the  astrometric  error  due  to  photon 
noise  is  negligible  and  the  major  errors  are  instrumental.  A 
detailed  discussion  of  one-  and  two-color  astrometry  will  be 
presented  in  separate  papers  (Colavita  et  al.,  1987;  Mozurkewich 
et  al.,  1988).  However,  the  astrometric  capability  of  the  instrument 
affects  all  other  modes  of  operation  of  the  instrument. 

With  the  known  coordinates  of  FK5  stars  it  is  possible  to 
derive  the  baseline  vector  and  delay  offset  by  measuring  the  fringe 
position  of  four  stars.  This  derived  baseline  is  as  accurate  as  the 
coordinates  of  the  FK5  stars  used.  An  accurate  baseline  vector  in 
turn  enables  the  instrument  to  find  the  fringe  faster  and  to  increase 
the  number  of  observations  from  several  per  night  to  several 
hundred  per  night.  This  astrometric  capability  is  the  key  to 
operation  of  the  instrument  as  an  absolute  interferometer,  de¬ 
scribed  next. 

3.2.  Absolute  interferometer 

Within  the  next  year  we  plan  to  install  the  software  for  real  time 
baseline  calibration.  In  this  mode,  the  interferometer  will  rapidly 
switch  between  4-6  FK5  stars  to  calibrate  the  baseline,  and  then 
point  blind  at  a  faint  astronomical  object  for  about  103s.  This 
cycle  will  be  repeated  to  yield  the  desired  integration  time  on  the 
faint  object.  Three  types  of  measurements  are  planned,  two  of 
them  in  the  near  future.  There  are  significant  advantages  to  phase 
coherent  operation  and  operation  in  the  absolute  mode  is  reserved 
for  objects  too  faint  for  the  fringe  tracker,  9  mag  for  10  cm  aper¬ 
tures.  All  three  types  of  measurements  planned  for  the  absolute 
interferometer  are  photon  starved  measurements.  Under  these 
conditions,  direct  phase  measurements  are  not  possible  and  the 
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Fig.  10.  Dispersed  fringe  detector  schematic 


The  signal-to-noise  ratio  of  the  above  visibility  estimator  has 
been  extensively  studied  both  for  amplitude  interferometry  and 
speckle  interferometry.  There  are  two  regions  of  interest,  depend¬ 
ing  on  photon  rate.  The  dependence  of  integration  time  needed  to 
achieve  a  given  signal-to-noise  ratio,  given  a  photon  flux  TV  per  r0 
aperture,  per  r0  integration  time,  per  spectral  channel  is  given  by 

Integration  time  oc  TV-1,  TV>  1,  TV >  Dark, 

ocTV-2,  1  >  TV >  Dark.  (11) 

For  stellar  diameter  measurements,  the  numerator  and  de¬ 
nominator  of  the  expression  in  Eq.  (10)  must  be  evaluated  every 
10  ms  for  each  spectral  channel,  typically  60  channels  of  50  A 
bandwidth  for  a  12-m  baseline.  If  collecting  apertures  greater  than 
r0  (10  cm)  are  used,  there  must  be  a  spectrometer  for  each  r0 
subaperture.  In  photon  starved  operation  there  is  no  hard 
magnitude  limit  as  there  is  in  active  phase  coherent  operation.  As  a 
practical  limit  we  require  a  signal-to-noise  ratio  of  30  in  104s  of 
integration.  Assuming  an  unresolved  star,  the  required  photon 
flux  is  1.6  10-2  photons  in  each  of  60  spectral  channels,  per  r0 
aperture  per  10  ms.  This  corresponds  to  a  magnitude  limit  of  14.3 
for  10-cm  apertures.  With  1-m  apertures,  the  limiting  magnitude  is 
about  16.8  (see  Appendix).  Accurate  visibility  measurements  will 
require  much  higher  signal-to-noise  ratios  and  the  magnitude  limit 
will  degrade.  The  optical  efficiency  of  the  spectrometer  was 
assumed  to  be  100%  and  the  quantum  efficiency  of  the  imaging 
photon  counter  was  assumed  to  be  10%.  For  these  reasons,  these 
estimates  should  be  considered  optimistic. 


only  quantities  that  can  be  measured  are  amplitude,  group  delay, 
and  closure  phase.  Photon  starved  conditions  occur  when  the 
average  number  of  photons  per  r0  aperture,  per  r0  integration 
time,  is  less  than  1  in  the  coherent  optical  bandpass. 

In  open  loop  operation,  a  12-m  interferometer  will  exhibit  rms 
phase  fluctuations  of  approximately  202  at  0.6  pm  for  1"  seeing. 
The  fringe  detector  must  be  able  to  accomodate  this  lack  of 
temporal  coherence.  On  the  other  hand,  to  maximize  sensitivity, 
light  from  all  wavelengths  should  be  used  for  interferometric 
measurements.  The  basic  fringe  detector  (Fig.  10)  consists  of  a 
prism  spectrometer  with  an  imaging  photon  counting  detector. 
The  required  resolution  of  the  spectrometer  is  of  the  order  of 

fr  \5/6 

,  (9) 

where  2  is  the  wavelength  of  observation,  L  the  length  of  the 
baseline,  and  r0oc26/5  is  equal  to  about  10cm  at  0.6pm  for  1" 
seeing.  Conservatively,  50  A  channels  are  assumed  in  the  calcu¬ 
lations  below.  In  the  following  sections  we  briefly  describe  the 
signal  processing  of  the  output  of  the  detector  for  the  three 
different  types  of  photon  starved  measurements. 


3.2.1.  Photon-starved  amplitude  measurements 


The  estimator  for  fringe  amplitude  must  be  modified  to  take  into 
account  a  bias  that  is  important  at  low  photon  rates.  An  unbiased 
estimator  for  visibility  squared  is  given  by 


7T2  r<(^-c)2  +  (^-D)2-A> 
2  [  < TV  —  Dark)2 


(10) 


where  N=A  +  B+C  +  D  =  total  photon  flux,  Dark  =  dark 
count  and  <  >  denotes  time  average. 


3.2.2.  Photon-starved  group  delay  astrometry 

The  signal  from  the  spectrometer  can  be  processed  in  a  different 
way  to  yield  astrometric  information  at  faint  magnitudes.  When 
the  optical  path  difference  in  the  two  arms  of  the  interferometer  is 
different  by  many  wavelengths,  say  3  pm,  there  will  be  construc¬ 
tive  interference  at  X  =  0.5  pm  and  0.6  pm,  and  destructive  inter¬ 
ference  in  between.  More  precisely,  the  spectrometer  output  will 
exhibit  a  fringe  pattern  due  to  the  interference  of  light  from  the 
two  siderostats  as  given  by 

/(v)  =  I0  (v)  (1  +  sin27z;  vt) ,  (12) 

where  I0  (v)  is  the  spectrum  of  the  star,  v  is  optical  frequency  (1/2) 
and  r  is  the  optical  path  difference  in  the  two  arms  of  the 
interferometer,  i.e.,  the  group  delay.  Equation  (12)  has  a  sine 
rather  than  a  cosine  because  the  Mark  III  uses  a  dielectric 
beamsplitter  to  combine  the  beams.  The  Fourier  transform  of  the 
“fringed”  spectrum  /(v)  will  have  a  delta  function  at  r.  At  this 
delay,  the  Fourier  transform  has  both  an  amplitude  and  phase.  In 
the  photon  starved  mode,  where 

§I(v)dv  <  1  photon/ro/10ms ,  (13) 

the  phase  is  not  measurable. 

The  basic  signal  processing  algorithm  is  to  average  the  power 
spectrum  of  /(v).  As  in  the  previous  section,  each  r0  subaperture 
must  be  separately  processed  if  the  collecting  apertures  are  large. 
Since  r  changes  due  to  turbulence,  a  time  average  of  the  power 
spectrum  is  a  gaussian  distribution  of  delta  functions.  For  a  12-m 
baseline,  the  standard  deviation  of  this  gaussian  is  ~12  pm  for 
1"  seeing.  A  SNR  of  30  in  this  case  means  that  the  peak  in  the 
gaussian  is  approximately  30  times  the  rms  fluctuations  in  the  rest 
of  the  power  spectrum.  For  10-cm  apertures,  the  magnitude  limit 
is  approximately  14.3  mag.  For  1-m  apertures,  the  corresponding 
limit  is  16.8  (see  Appendix). 
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3.2.3.  Photon-starved  closure  phase 


!  Closure  phase  measurements  were  first  made  with  radio  inter- 
[  ferometers  (Jennison,  1958),  and  subsequently  suggested  for 
[  optical  interferometers  (Rogstad,  1968).  Recently,  closure  phase 
>  measurements  were  made  at  visible  wavelengths  by  masking  the 
!  aperture  of  a  large  telescope  (Baldwin,  1986;  Readhead,  1988)  at 
1  photon  rich  levels.  In  closure  phase  measurements,  three  inter¬ 
ferometers,  whose  baselines  form  a  closed  triangle,  measure  three 
complex  visibilities  corrupted  by  turbulence.  However,  the  closure 
phase,  the  sum  of  the  phases  of  the  three  complex  visibilities,  is 
independent  of  the  atmosphere.  In  photon  starved  measurements, 
raw  phases  cannot  be  measured.  Instead,  the  product  of  the 
complex  visibilities  is  averaged,  i.e., 


triple  product  =  (At  ei<Pi  A2ei(p2  A3< 


3>, 


(14) 


where  A±  ei<Pi  is  the  complex  visibility  from  interferometer  1.  The 
importance  of  averaging  the  complex  visibilities  comes  from  the 
fact  that  the  product  of  complex  visibilities  is  not  biased  in  the 
photon  starved  region.  The  closure  phase  is  simply  the  imaginary 
part  of  the  log  of  the  triple  product  (see  Appendix).  The  triple 
product  was  first  used  in  high  resolution  optical  imaging  for 
reconstruction  of  speckle  images  (Weigelt,  1977),  a  technique 
called  speckle  masking.  This  technique  extends  the  sensitivity  of 
closure  phase  measurements  more  than  an  order  of  magnitude 
beyond  conventional  photon  rich  measurements. 

For  long  baseline  interferometers,  the  triple  product  must  be 
calculated  every  10ms,  for  each  spectral  channel,  and  each  rQ 
aperture.  The  signal-to-noise  ratio  (SNR)  of  the  triple  product  is 
qualitatively  different  from  the  SNR  for  amplitude  measure¬ 
ments.  At  high  photon  fluxes,  A>  1,  the  SNR  of  all  these 
quantities  varies  as  j /a.  At  low  photon  fluxes,  A  <  1 ,  the  SNR  for 
amplitude  measurements  varies  as  A,  and  that  of  the  triple 
product  varies  as  A3/2  (see  Appendix).  Again,  assuming  an  SNR 
of  30  in  104  s  (A<j>  measured  to  2/1 80),  the  magnitude  limit  with  10- 
cm  apertures  is  approximately  12.5  mag,  or  a  photon  flux  of  0.08 
photons  per  r0,  per  10  ms,  per  50  A  bandwidth.  The  limiting 
magnitude  for  1-m  apertures  is  14.2. 

Although  three  siderostats  will  be  in  operation  in  1987,  closure 
phase  measurements  are  not  planned  for  the  immediate  future.  In 
order  to  operate  three  interferometers  simultaneously  the  second 
delay  line  must  be  converted  to  high  precision  operation.  Oper¬ 
ation  as  an  absolute  interferometer  will  then  enable  closure  phase 
measurements  to  be  made  with  a  one  to  two  order-of-magnitude 
increase  in  sensitivity  over  active  phase  coherent  operation. 


4.  Conclusion 

In  the  past  seventy  years  a  number  of  stellar  interferometers  have 
been  proposed  and  built.  Starting  with  Michelson,  the  instruments 
have  gotten  progressively  more  complex  and  sophisticated 
(Koechlinetal.,  1985;  Labeyrieet  al.,  1986;  Davis  etal.,  1986).  The 
Mark  III  incorporates  many  of  the  ideas  proven  by  these  previous 
instruments.  The  major  differences  between  this  interferometer 
and  others  reside  primarily  in  its  operational  characteristics  -  it  is 
fully  automated,  it  is  easy  to  operate,  and  it  is  mechanically 
accurate  and  thermally  stable.  It  is  also  capable  of  making  literally 
hundreds  of  observations  in  one  night;  both  astrometric  and 
stellar  diameter  measurements  are  possible.  Last  of  all,  this 
instrument  has  the  potential  to  operate  as  an  absolute  inter¬ 
ferometer  which,  from  an  astronomical  point  of  view,  holds 
tremendous  promise  for  observation  of  faint  objects. 


To  fully  exploit  the  possibilities  of  absolute  interferometry, 
and  to  construct  images  of  faint  extragalactic  objects  at  milliarcsec 
resolution,  even  more  complex  interferometers  with  large,  mul¬ 
tiple  apertures  will  be  needed.  During  the  next  several  years  the 
Mark  III  should  demonstrate  much  of  the  necessary  technology, 
as  well  as  produce  some  of  the  scientific  results  that  would 
motivate  the  construction  of  a  large  optical  array. 
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Appendix:  photon-starved  interferometry 

The  photon-starved  signal-to-noise  ratio  for  the  three  inter¬ 
ferometric  quantities:  amplitude,  group  delay,  and  closure  phase 
are  discussed  in  this  Appendix.  We  present  limiting  magnitudes 
for  these  quantities  extrapolated  from  our  experience  with  the 
Mark  III  interferometer.  Recently,  there  has  been  considerable 
interest  in  the  use  of  light  scattered  from  the  ionosphere  from  a 
pulsed  laser  in  order  to  provide  a  probe  of  atmospheric  turbulence 
to  operate  an  active-optics  system  (Foy  and  Labeyrie,  1986).  The 
application  of  this  technique  is  not  considered  in  our  calculations. 

Estimates  of  amplitude,  group  delay,  and  closure  phase  each 
have  slightly  different  magnitude  limits.  Qualitatively,  different 
effects  degrade  the  sensitivity  for  the  different  quantities.  In  all 
cases,  we  assume  that  the  fringe  detector  is  a  dispersive  spec¬ 
trometer  followed  by  an  imaging  photon-counting  detector.  We 
also  assume  that  the  resolution  of  the  spectrometer  is  a  few  times 
greater  than  what  is  needed  to  maintain  high  fringe  contrast  in  the 
presence  of  the  optical  path  fluctuations  caused  by  atmospheric 
turbulence.  Magnitude  limits  are  calculated  assuming  an  unre¬ 
solved  star,  a  12-m  baseline,  10-cm  collecting  apertures,  a  10  ms 
coherent  integration  time,  a  60-channel,  50  A  resolution  spec¬ 
trometer,  unit  fringe  visibility,  a  total  integration  time  of  104  s,  and 
a  final  signal-to-noise  ratio  of  30.  The  detector  dark  count  is 
neglected,  and  a  total  system  efficiency  of  0.1  is  assumed.  We 
ignore  the  reduction  in  fringe  visibility  and  hence  signal-to-noise 
ratio  due  to  the  use  of  coherent  apertures  of  10  cm  and  integration 
times  of  10  ms.  In  practice,  the  seeing  would  have  to  be  somewhat 
better  than  1"  for  these  magnitude  limits  to  be  achieved.  Scaling 
formulae  are  given  for  the  case  of  different  seeing  conditions. 

At  high  photon  fluxes,  when  the  number  of  photons  per  r0 
aperture,  per  atmospheric  coherence  time  t0  ,  per  spectral  channel 
AX,  is  much  greater  than  one,  the  signal-to-noise  ratio  is  propor¬ 
tional  to  the  square  root  of  the  photon  flux  for  estimates  of  all 
three  interferometric  quantities.  We  define  the  quantities  A  and  M 
as 

A  =  No.  of  detected  photons  per  frame,  i.e.,  per  r0,  t0,  AX, 
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where  D  =  diameter  of  telescope,  T  =  total  integration  time,  and 
BW  =  total  optical  bandwidth  of  the  spectrometer  (assumed  0.4  to 
0.7  pm).  The  terms  photon-rich  and  photon-starved  apply  to  the 
cases  A  >  1  and  A  <  1 ,  respectively. 

The  signal-to-noise  ratio  for  fringe  amplitude  measurements 
has  been  extensively  studied  for  many  types  of  interferometers.  It 
is  usually  more  convenient  to  consider  estimates  of  V 2  rather 
than  of  V.  The  signal-to-noise  ratio  for  estimates  of  V2  can  be 
written  (Tango  and  Twiss,  1980) 


SNR(V2)  =  ^M112  NV2 


1  +  ^NV2 


-1/2 


(Al) 


For  the  case  of  a  fringe-scanning  interferometer  with  4  bins  per 
scan,  A  is  replaced  by  (8/7t2)  A.  The  qualitative  feature  of  this 
expression  is  that  at  low  photon  fluxes,  A  <  1 ,  the  signal-to-noise 
ratio  is  proportional  to  A,  rather  than  to  A 1/2 .  Thus,  if  a  star  is 
dimmer  by  a  factor  of  10,  the  number  of  frames  M  must  increase 
by  a  factor  of  100  to  give  the  same  signal-to-noise  ratio. 

The  noise  behavior  of  the  group-delay  estimate  is  similar  to  the 
V2  estimate.  The  “fringed”  spectrum  as  detected  by  the  spec¬ 
trometer  is 


/(v)  =  /0(v)(l  +sin27rvT'),  (A2) 

where  I0(v)  is  the  spectrum  of  the  star,  v  is  optical  frequency,  and  x' 
is  the  group  delay.  The  group  delay  changes  due  to  atmospheric 
turbulence.  Over  long  time  intervals,  x'  is  a  gaussian-distributed 
random  variable  with  mean  ,  determined  by  the  position  of  the 
star  and  optical  delay  line,  and  variance  a2 ,  which  depends  on  the 
seeing  and  baseline  length.  The  variance  is  given  by 

a2  =  6.88  (Xjln)2  (L/r0)513  ,  (A3) 


where  L  is  the  baseline  length  and  r0  is  the  atmospheric  coherence 
diameter. 

Each  coherence  time  r0,  the  power  spectrum,  P( r),  of  /(v)  is 
calculated.  The  power  spectrum  is  then  averaged  for  the  length  of 
the  observation.  There  are  three  important  features:  a  peak  at  the 
origin  attributable  to  the  Fourier  transform  of  70  (v),  a  uniform 
bias  due  to  photon  noise,  and  a  gaussian  distribution  at  the  mean 
group  delay  xl  due  to  atmospheric  turbulence.  The  averaged 
power  spectrum  is  given  by 

<P(t)>  =  <| \e2niVT I0(v)  (1  +sin2^vT')  dv |2>,  (A 4) 

where  <  )  represents  an  ensemble  average  over  both  group  delay 
x'  and  photon  statistics. 

A  better  estimator  is  the  unbiased,  normalized  power  spectrum 
Q  (t),  given  by 


GW  = 


<P(T)-fl(y)dv> 

<J/(v)t/v>2 


(A5) 


where  Ns  =  j  I(v)dv  is  the  total  photon  flux  in  a  10  ms  interval. 
Comparison  of  this  expression  with  the  estimator  (11)  shows  that 
(A  5)  is  just  the  V2  estimator  evaluated  at  a  particular  group  delay. 
Thus,  the  performance  of  the  group  delay  estimator  can  be 
calculated  using  (Al). 

Given  a  spectrometer  with  a  bandwidth  of  0.4  pm  to  0.7  pm, 
the  delay  resolution  is  Ax  —  (1/0.4  — 1/0.7) -1  ~  1  pm.  Since  the 
group  delay  fluctuates  with  a  standard  deviation  of  erT,  the 
gaussian  distribution  of  group  delay  in  Q  (r)  will  be  spread  over 
approximately  aJAx  bins,  while  the  amplitude  of  the  gaussian  will 
be  reduced  by  approximately  crJAx.  Photon-starved  statistics 
apply  for  (1/2)  Ns  V2  <  aJAx,  and  using  (A  1),  the  signal-to-noise 


ratio  for  an  estimate  of  the  amplitude  of  the  gaussian  using  a 
matched  filter  is  given  by 

N  V2M112 

SNR(Q (t))  = 

where  As  is  the  total  photon  flux  in  the  spectrum.  In  (A  6)  sampling 
losses  are  neglected,  and  it  is  assumed  that  ox^>  Ax.  Using  (9)  and 
(A3)  to  write  (A 6)  in  terms  of  A ,  rather  than  Ns ,  yields  a  signal-to- 
noise  ratio  for  group  delay  estimates  which  is  nearly  identical  to 
the  signal-to-noise  ratio  of  a  conventional  amplitude  estimate. 

For  the  case  of  group-delay,  an  estimate  of  the  centroid  of  the 
group-delay  distribution  is  desired.  The  standard  deviation  of 
such  an  estimate  can  be  written 


where  the  SNR  is  given  by  (A  6).  A  Monte  Carlo  simulation  was 
performed  which  verified  the  above  analytically  derived 
expressions. 

Under  photon-starved  conditions,  the  measurement  of  closure 
phase  is  qualitatively  different  from  the  measurement  of  ampli¬ 
tude.  The  definition  of  the  triple  product,  Eq.  (14),  is  repeated 
below: 

triple  product  =  ei<Pl  A2  ei<p2  A3  ei(p 3>  .  (A 8) 

In  photon-counting  measurements,  if  zero  photons  are  detected  in 
any  one  of  the  three  interferometers,  the  triple  product  is  zero.  For 
a  frame  to  contain  some  information  on  the  closure  phase,  at  least 
one  photon  must  be  detected  in  all  three  interferometers.  The 
Poisson  probability  function  is  given  by 

Nme~N 

Pim)=^T’  (A9) 


where  P(m)  is  the  probability  of  detecting  m  photons  given  an 
average  arrival  rate  of  A.  For  N  <|  1  the  probability  of  detecting 
one  photon  is  N.  Thus,  the  probability  of  detecting  one  photon  in 
all  three  interferometers  is  A3,  so  that  the  variance  for  estimates  of 
closure  phase  should  be  proportional  to  A-3  under  photon 
starved  conditions. 

More  precisely,  the  variance  of  the  closure  phase  calculated  via 
the  triple  product  can  be  written 


for  a  signal-to-noise  ratio  SNR  =  l/o^ .  In  this  expression  A  is  the 
average  number  of  detected  photons  in  each  interferometer.  Equal 
visibilities  V  have  been  assumed  for  each  of  the  three  inter¬ 
ferometers.  Thus,  at  high  photon  fluxes,  the  variance  is  three  times 
the  variance  in  a  photon-rich  phase  measurement.  At  low  photon 
fluxes,  the  variance  is  proportional  to  A-3,  as  given  by  the  last 
term  of  (A  10).  For  this  term,  V  is  the  geometric  mean  of  the 
visibilities  in  the  three  interferometers.  A  Monte  Carlo  simulation 
for  closure  phase  measurements  was  also  conducted,  which 
verified  the  above  expression. 

The  integration  time  needed  to  achieve  a  signal-to-noise  ratio 
of  30  versus  photon  flux  A  is  given  in  Fig.  1 1  for  measurements  of 
visibility  squared  and  group  delay,  and  closure  phase.  For 
astronomical  measurements,  a  photon  flux  corresponding  to  an 
SNR  of  30  on  an  unresolved  star  allows  a  diameter  measurement 
of  4%  on  a  partially-resolved  source.  An  SNR  of  30  for  group 
delay  corresponds  to  an  astrometric  measurement  of  0.6  pm  in 
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delay,  or  0"01  in  angle.  An  SNR  of  30  for  the  triple  product 
corresponds  to  an  error  of  2/180  in  closure  phase.  The  limiting 
magnitudes  for  a  single-r0  system  are  approximately  14.3  and  12.5 
for  visibility/group  delay  and  closure  phase.  With  1-m  aterpures, 
the  limiting  magnitudes  are  approximately  16.8  and  14.2, 
respectively. 

If  r0  and  r0  are  different  from  the  nominal  values  assumed,  the 
magnitude  limits  will  change.  In  general,  an  instrument  will  have  a 
fixed  collecting  aperture.  Thus,  if  r0  increases,  the  number  of 
photons,  N ,  per  r0  increases,  but  the  number  of  frames,  M, 
decreases.  Similarly,  given  a  fixed  104s  for  a  measurement,  an 
increase  in  t0  will  also  increase  N  and  decrease  M .  If  the  signal-to- 
noise  ratio  were  dependent  on  N  and  M  evenly,  the  final  SNR 
would  be  independent  of  seeing.  This  is  the  case  for  photon-rich 
measurements  where  the  SNR  is  proportional  to  j/AM . 

For  photon-starved  measurements,  this  is  not  the  case.  The 
signal-to-noise  ratio  for  amplitude  and  delay  measurements 
follow  the  scaling  relation  AM1/2,  while  the  SNR  for  closure  phase 
follows  the  relation  N312  M1/2.  Two  additional  facts  should  be 
considered.  As  r0  improves,  the  rms  path  length  fluctuation 
between  the  apertures  decreases  proportionally.  Hence,  the  num¬ 
ber  of  required  spectral  channels  that  must  be  separately  analyzed 
decreases  according  to  (9).  In  addition,  since  the  temporal 
fluctuations  of  phase  are  usually  modeled  by  assuming  frozen 
turbulence  translating  past  the  instrument  at  a  constant  wind 
speed,  Tq  should  increase  with  r0 .  Thus,  the  scaling  of  N  and  M 
with  r0  is  approximately  given  by 

Aocr£,  (All) 

Macro4.  (A12) 


The  scaling  of  the  magnitude  limits  for  photon-starved  measure¬ 
ments  with  seeing  is  then  given  by 

Amplitude/delay  sensittivity  oc  ,  (A  1 3) 

Phase  closure  sensitivity  oc  r®/3  .  (A  14) 

These  results  should  be  compared  to  the  case  of  a  wideband  phase 
tracking  system.  In  this  case,  only  N ,  not  Mis  important,  and  since 
the  fringe  is  tracked,  there  is  no  need  for  the  spectrometer.  Thus, 
the  sensitivity  scales  with  seeing  as  rl . 
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